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Abstract 

A second quantised theory of electrons and positrons in a deep time-dependent 
potential well is discussed. It is shown that positron production from the well 
is a natural consequence of Dirac's hole theory when the strength of the well 
becomes supercritical. A formalism is developed whereby the amplitude for 
emission of a positron of a given momentum can be calculated. The difference 
between positron production and electron-positron pair production is demon- 
strated. Considerations of the vacuum charge and of Levinson's theorem are 
required for a full description of the problem. 
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1. Introduction 



Gershtein and Zeldovich [1] showed in 1969 that spontaneous positron production was 
possible when two heavy bare nuclei with total charge Z greater than some critical 
value Z c collided with each other. Z c is the value of Z for which the lowest IS state 
of the hydrogenlike atom with total charge Z distributed on a nucleus of radius R has 
energy E = —m, where m is the electron mass. Pomeranchuk and Smorodinskii [2] 
calculated Z c = 200 for a nucleus assuming the electrostatic potential is constant for 
r<R, while Greiner and Pieper [3] calculated Z c = 169 for a homogeneous spherical 
charge distribution. 

Zeldovich and Popov [4] subsequently reviewed the theoretical basis of this process. 
(For a recent overview of the problem and of other related topics the reader may find 
the monograph [5] useful.) They were able to set up a self-consistent picture of spon- 
taneous positron production but from the viewpoint of a first-quantised Dirac theory. 
They state, however, that "at Z > Z c we have in principle a many-body problem. 
This constitutes the third stage of the investigation of the Z > Z c situation, and has 
not been completed as yet. In an exact formulation, it is necessary to consider the 
equations of the electron-positron wave field with creation and annihilation operators, 
and to use the second quantisation formalism". They go on to say that they predict 
the results of experiments on the basis of single particle solutions of the Dirac equa- 
tion but that "such conclusions and predictions require a justification, which can be 
obtained only from the exact many-body theory" 

This is the problem we have set ourselves in this paper for as far as we know it 
has not yet been tackled. For simplicity we have considered spontaneous fermionic 
production in the simplest supercritical external field we could think of, which is a 
sufficiently deep one-dimensional square well. We shall find that many of Zeldovich 
and Popov's results can indeed be justified. We shall also see that the preparation of 
the supercritical potential is an essential feature in the problem and that while it is 
true that a static supercritical square well does not radiate [6], it is also uninteresting 
because in any experimental situation it is necessary to create the supercritical poten- 
tial from sub-critical potentials (by, for example, colliding two heavy ions together) 
and thus the time-dependence of the external field is an essential ingredient. 

We first have to elucidate various questions relating to the Dirac equation in one 
dimension. In section 2 we consider a Dirac particle in a square well, write down 
the spectrum of bound states and give expressions for the scattering amplitude and 
the phase shift. We also present the conditions for criticality. In section 3 we give 
the second quantized version of the theory paying particular attention to possible 
conventions related to the labelling of states. We review the concept of the vacuum 
charge that plays a central role in this discussion. In section 4 we derive Levinson's 
theorem for the Dirac equation in one dimension and express the vacuum charge in 
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terms of phase shifts. In section 5 we illustrate some of the above concepts in the 
case of a 5 function potential. In section 6 we apply the methods and results of the 
previous sections to the particular case of an external field which makes the transition 
from slightly subcritical to slightly supercritical to slightly subcritical and show in 
detail that fermionic radiation can occur. For this it is necessary to arrange for the 
ground state of the well to be initially empty, but we demonstrate how we end up 
with a filled bound state and a free fermion (which we identify for convenience with a 
positron). Second quantization and use of the Bogoliubov transformation are essential 
in this treatment. 

2. The One Particle Dirac Equation in One Dimension 
2.1. The Free Dirac Equation. 

We adopt the convention 7o = 0-3, 7i = icri- The above choice agrees with 7i7j + 7j7i = 
2gij. The free Dirac Hamiltonian in one dimension is then 

H = -cr 2 p + cr 3 m 
and the Dirac equation takes the form 

(a 1 ^-Ea 3 + m)V = (2.1) 

In what follows k stands for the wavevector, k for its magnitude and e = \E\ = 
+\ / k 2 + m 2 . We try a plane wave of the form 

\ elkx _ iEt ^ ^ 

and substitute in (|2.1|) . The equation is satisfied by A = ik, B = E — m where 
E = ±£. Note that in one dimension the negative energy solution is obtainable from 
the positive energy one simply by replacing e by — e unlike what happens in three 
dimensions. We normalize the 2 dimensional spinor in ( |2.2| ) by requiring / dxu^u = 1. 
The positive energy solutions have the form 

and the negative energy ones 

NJe) ( lk 

\ — £ — m 



^ikx+iet ^ 4^) 
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where N±(e) are normalization factors. In what follows we will either consider the 
particle in a box of length 2L and imply periodic boundary conditions at x = —L and 
x = L or use continuum normalization. In the former case 

N + (e) = — , L , M( e ) ' 



V2L^/2e{e - m) V2L^2e{e + m) 

and in the latter 

Nje) = _ , 1 , i\L( £ ) = — 



y/2iryj2e{e - m) ^/2-n^2e[e + m) 

Finally we quote the transformation of a wavefunction under parity: 

ijj'(x,t)=a^(-x,t) (2.5) 

2.2. The Dirac Particle in a Square Well: Bound States. 

Consider a square well of depth —V extending from x = —a to x = a. The electron 
Hamiltonian is 

Hq = —&2P — V + &3 m 

where V is the zeroth component of A^. The electron charge is taken to be -1. The 
Dirac equation reads 

(a l ^--(E + V)a 3 + m)ij = (2.6) 

Define 



p= y(E+ V) 2 — m 2 (2.7) 
Inside the well try the wavefunction 

/ Acospx \ 
I B sin px J 

According to (|2.5| ) this describes a state even under parity. Substituting in ( |2~T| ) 
we determine A, B (modulo a factor): 

_/ pcospx \ 

Weven ^ ( E + V - m) Sin P X J 1 ' 

for —a < x < a. For x > a it can easily be checked that 

= 8 ( E ~^ m ) e— (2.9) 

with 

k = V-E 2 - m 2 (2.10) 
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satisfies the Dirac equation. Matching ( |2.8| ) and ( |2.9| ) at x = a and eliminating s we 
obtain the equation that determines the spectrum of the even bound states: 



tanap 



(m- E)(E + V + m) 
\l (m + E)(E + V -m) 



(2.11) 



Note that it is implicit in (2.7) that for a subcritical potential (i.e. V < 2m; see 
below) permissible bound states (even or odd) satisfy E > m — V . In a similar way 
we can determine the wavefunctions of the odd bound states. It turns out that for 
—a < x < a 

, _ / — psinpx \ 
%dd - y (E + V -m)cospx ) 

and for x > a 

.1. J I ^ 1 „ — KX 



(2.12) 



E — m 



(2.13) 



Matching ( 2.12j ) and (p. 13 ) at x = a and eliminating s' we obtain the equation that 
determines the spectrum of the odd bound states: 



tanap 



{m + E)(E + V-m) 



\ (m- E)(E + V + m) 
These results are well known and can be found, for example, in reference [3]. 



(2.14) 



2.3. The Dirac Particle in a Square Well: Scattering States. 

Consider a wave incident from the left. The corresponding wavefunction is 



for x < a, 



ik 
E — m 



ip 



ikx 



B 



—ik 
E — m 



-ikx 



\ E + V — m 



—ip 
E + V - m 



-ipx 



for —a < x < a and 



ik 
E — m 



Akx 



(2.15) 



(2.16) 



(2.17) 



for x > a. To calculate the wavefunction we require continuity at x = —a and at 
x = a. The coefficients B,C,D,F are calculated in 0. We only give the results for 
the phase shift. Define 

F(E) = \F(E)\e i5{E) 



and 



kE + V + m 



7 = - 



p E + m 
5 



(2.18) 



Then 

6(E) = arctan tan2poj - 2ka (2.19) 

In what follows the phase shift will sometimes be written with two arguments, i.e. 
S(E,V), the second referring to the particular value of the potential. Whenever we 
write 5 (E, 0) we understand the limit 6 (E, V — > 0). The algebra leading to (2.19) is 
right regardless of the sign of the energy. Recalling that e = \E\ the phase shifts for 
positive and negative energy will be written as 5 + (e, V) , <5_ (e, V) respectively. The 
phase shift is defined modulo a multiple of 7T. We choose 

5(E,0) = (2.20) 

(if V = then k ~ p , 7 — > 1 and (|2.19| ) satisfies (2.19)). We also write 



5(±oo, V) = ±2Va (2.21) 

Indeed if E +00 (with V, m fixed) then k ~ E,p ~ E, ( |2.18| ) yields 7 ~ 1 and from 
( |2.19| ) 5 ~ 2pa — 2A;a — > 2Va. Similarly for — > —00 , 7 ~ 1 and 5 ~ 2pa — 2fca 



2V^a. Note that (|2.21|) satisfies 5(oo) + 5(— 00) = as shown in 3 dimensions by Ma 



and Ni [8]). We also have 

/oo 
V(x)dx 
-00 

which is analagous to the 3-dimensional result 

roo 

6 (±00) = =F / V(r)dr 
Jo 

It must be emphasized that this result is in contrast to what happens in nonrelativistic 
quantum mechanics where the phase shift at high energy always tends to zero (modulo 
a multiple of it). The absence of multiples of n in the right hand side of (2.20), (2.21) 
is not a matter of convention. That would be the case had we specified e.g. 5 (+00, V) 
only. Having however fixed the phase shift at one end of the energy range the rest 
should follow. It is better to regard (2.20), (2.21) as dictated by physical considerations 
as well. For V — » we expect the results of the Dirac equation to coincide with those 
of the Schroedinger equation and (2.20), (2.21) satisfy this requirement (without any 
multiples of 71"). Notice that the negative energy phase shift makes sense in the context 
of nonrelativistic quantum mechanics as well: it corresponds to the scattering of a 
positron of momentum — k. 

We now wish to make a statement concerning the threshold values <5(±m, V). For 
a certain fixed V consider positive ( negative) energy states very near threshold, i.e. 
E — > m+ (E — > —m—). Then the quantity 7 defined in fl2.18|) has small absolute 



value and is positive (negative). Hence from ( 2.19| ) we conclude that 



S(m, V) = ^ + n(V)7i (2.22) 
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5{-m,V) = ~ + n , (V)ir (2.23) 

where n (V) , n' (V) are integers. Notice that at this stage no statement is made 
about n (V) , n' (V). In particular the values at threshold n (0) , n' (0) are unknown; 
see section 4 for details. Although the above equations have the look of Levinson's 
theorem in nonrelativistic quantum mechanics (see e.g. [[]], p. 354) in the case of the 
Dirac equation it is the sum n + n! that is connected to the number of bound states 

7T 

[8]. We consider this point in more detail in section 4. However the occurrence of — 

in ( |2.22|) and (|2.23|) has the same origin as the — in the corresponding equation in |J, 

namely the presence of a half-bound state in one-dimension [9]. 

We are particularly interested in the case where there is no reflection and so the 
transmission coefficient is a maximum. These we call, following Bohm[10], transmis- 
sion resonances. Setting B = and matching the above expressions at the boundaries 
we get four equations and eliminating C, D, F we obtain the equation determining the 
spectrum of the resonances 

2pa = Ntc (2.24) 

This relation connecting wavelength and dimension of the well is familiar from the 
Fabry-Perot etalon. The similarity in the mathematical treatment of Fabry-Perot 
and a Dirac square well was first pointed out in reference [6]. We should emphasize, 
however, that transmission resonances are not resonances in the normal sense; the 
phase shift does not increase through n/2 as it does with a proper resonance. 

A further point concerning the relation between bound states and transmission 
resonances may be in order. Consider a transmission resonance of energy E with k 
the incident wavevector and p the wavevector between the walls. Using the boundary 
conditions one can easily show that 

^ = e i2p ' a (2.25) 



The above equation together with ( |2.24| ) yield C = ±D. Substituting in (|2.16|) we 



get that the transmission resonance wavefunction between the walls for N — 1 has 
the form ( |2.8| ). In other words if we consider an almost critical potential then the 
wavefunction of the lowest lying (even) bound state and the wavefunction between 
the walls of the upper (N = 1) resonance approach each other as the potential varies 
and the energies Eb oun d and E res approach — m from above and below respectively. 
Hence the overlap integral betwen the two wavefunctions is sizeable. 

We define positive energy spinor wavefunctions «( + )(k, x) given by Q2.15D , ( |2.16|) , 
fl2.17p multiplied by N + (e) and negative energy ones W(_)(k, x) replacing E by — E 



multiplied by N_(e). Note that u^u^ = ut )«(+) = and 

J dxu^(k, x)w( + )(k', x) = J dxu^(k, x)-U(_)(k', x) = 5(k — k') (2.26) 
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u^(k,x) = u(— k, x) for both kinds of spinors. These wavefunctions refer to scattering 
states. We also have to include bound states. The spinor wavefunction corresponding 
to the jth bound state is denoted by Uj(x) and can be read off Q2.8|), (|2.9|), ( |2.12| ), ( |2.13j ) 
(at present the notation does not distinguish between even and odd bound states). 
Uj is assumed to be normalized / dxUjUj = 1 and is orthogonal to the continuum 
wavefunctions. 

Since parity is a good quantum number it is convenient to introduce spinor wave- 
functions that transform under parity in a definite way. Using (|2.5|) define parity even 
and odd spinors corresponding to positive energies: 

u (+)e {k, x ) = ^ ( M (+)( k > x ) - 03«(+)(-k, x)) (2.27) 

u {+)o (k, x) = -j= (u {+) (k, x) + cr 3 w (+ ) (-k, xfj (2.28) 

Replacing the subscript (+) by (-) we obtain the corresponding expressions for negative 
energy spinors. 



2.4. Criticality 

Let us now address the notion of criticality that plays a crucial role in the discussion. 
It is quite clear that for a small value of the potential equation ( |2.11|) always has a 
solution corresponding to an even bound state. Suppose that the potential deepens 
gradually. At some stage the energy of the bound state crosses zero; the implications 
of this effect on the vacuum charge will be discussed later. For greater values of 
the potential the energy approaches — m and for a critical value Vf (the subscript 1 
refers to the fact that this is the first bound state that disappears) the bound state 
merges with the negative energy continuum states. Whether other bound states have 
appeared in the meantime is considered presently. The wavevector p between the walls 
is related to the potential by ( |2.7| ) and solving for V we get 



V = yV + m? - E (2.29) 

It is clear from (|2.11 ) that when 



ap = | (2.30) 

the even bound state is at — m and subsequently disappears. This corresponds to a 
critical value of the potential 



V^ = \l— + m? + m (2.31) 



Comparing Q2.30 ) with ( p. 24 ) we see that the bound state wavefunction goes over to 



the N = 1 transmission resonance wavefunction. On the other hand it is clear from 



S 



71 

(2.14) that when ap = — a new odd bound state appears at E = m. This corresponds 
to the value of the potential 



7T 2 



K^i = y^ + m 2 -m (2.32) 



When pa = it the second even bound state appears at 



/tt 2 

V eV en,2 = \l ~^ + m 2 - 111 (2.33) 

Again for pa = tt the first odd bound state disappears corresponding to the poten- 
tial 

and to the N = 2 transmission resonance. It should by now be clear how the pattern 
repeats itself. 
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3. Second Quantization. 

Let N + ,N_ be the number of positive and negative bound states respectively. We 
expand ip in terms of the continuous spectrum wavefunctions ( P-27 ), (|2.28| ) and of the 
bound states: 

ijj(x, t) = J2ki a e(k, 0)u( + ) e (k, x)e~ iEt + a a (k, 0)u {+)o (k, x)e~ iEt + 

+4{k, 0)« ( _ )e (fc, x)e iEt + cl{k, 0)u { - )o {k, x)e iEt }+ (3.1) 

+{Ef=i &,(0)^(x)e-^* + Z?=i 4{0) Uj {x)e- iE ^} 

Operators a), a create and annihilate travelling electrons; c\c are the corresponding 
ones for positrons. Operators bj (fcj) annihilate (create) bound electrons whereas dj 
(dj) annihilate (create) bound positrons. The use of the (f) in ( |3.1| ) is dictated by the 
sign of the exponential and conforms to current literature. The Hermitean conjugate 
expansion is 

V'tfo t) = T, k {at(k, 0)u\ +)e (k, x)e iEt + al(k, 0)u\ +)o (k, x)e iEt + 

+c e (k, 0)u\_ )e (k, x)e~ lEt + c {k, 0)u\_ )o (k, x)e-* Et }+ (3.2) 
+{Eft b](0)u](x)e iE ^ + Ef=i ^(OKOr)^} 

where we took into account the reality of Uj. The standard anticommutation relations 
are obeyed 

{a e>0 (k,t),al (k',t} = 5 kk '5 eo , {bi(t),b}(t)} =8^ 

(3.3) 

{c e>0 (k,t),cl (k',t} = 5 kk ,5 eo , {di(t),dj(t)} = 5^ 

We work in the Heisenberg picture throughout: The time dependence is carried by 
operators whereas state vectors are time independent. However basis ket vectors (and 
in particular the vacuum) are time dependent (see e.g. |TT|]). The vacuum |0) is 
defined by 

a e ,o(k) |0> = c efi {k) |0> = h |0> = di |0> = (3.4) 

It must be stressed however that the definition of the vacuum is to some extent a 
matter of convention. A complete statement on the physics of the problem is made 
when one defines both the vacuum and also the state vector of the system. We will 
elaborate this point later on. 
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The total charge is defined by (N.B. according to our conventions the electron 
charge is —1) 

Q(t) = / dxp(x,t) = — / dx ip\x, t), ip(x, t) (3.5) 



2 

Substituting ( J3.1|) , ( |3.2| ) and using ( ^.3[ ) we get 



Q = Q norm + Qo (3.6) 



where 



and 



Qnorm = Efc{-4(&; *0 a e(&, *) + 4(fc, t)c e (k, t) - dl(k, t)a (k, t)~ 

+cl(k, t)c (k, t) - T!j=i bjbj + Ey=i djdj} 

Qo = - < ^(states with J5 > 0) - ^(states with E < 0) > 
^ I k k ) 



(3.7) 



(3i 



Expressions (|3.7p , ( |3.8| ) are rather formal. To be precise we have to calculate phase 
shifts and density of states and transform the sums to integrals. This is done in section 
4 and illustrated in section 5. Given the definition ( |3.4| ) of the vacuum we immediately 

get 

(0|Q|0) = Q (3.9) 

i.e. the vacuum charge turns out to be the spectral asymmetry of the Hamiltonian. 
Note that we would not have obtained the connection between Qo and the vacuum 
asymmetry had we not identified bound states with E < as positrons or not used the 
commutator in ( |3.5| ). Note also that Q clearly depends on the potential. The same 
applies to Q norm although the notation does not indicate this. The wavefunctions 
that appear in (3.1), (3.2) depend on the strength of the potential. By implication 
the same is true for the creation and annihilation operators and thus for Q nor m- 

To see precisely how the above formalism should be used let us consider the electron 
field interacting with a time dependent potential V(x,t). The Hamiltonian is given 
by 

If r i 

H = H(freeDirac) - - / dxV(x, t) ip^x, t), ^(x, t) (3.10) 



2 

It is quite clear that (|3.10| ) conserves Q. In order to avoid extraneous complications 
let us assume for the moment that the potential is relatively weak and that there is 
only one bound state at energy Ej, slightly above zero , i.e. N + = 1, i\L = 0. Denote 
by Qo,init the initial vacuum charge. In the present case ( |3.1| ) reads 

ij)(x, t) = J2k{ a e(k, 0)u {+)e (k, x)e' iEt + a D (k, 0)u (+)o (k, x)e~ iEt + 

(3.11) 

+4(k, 0)u ( _ )e (fc, x)e lEt + cl(k, 0)u Ho (k, x)e lEt } + b(0)u(x)e- lEbt 
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Before going any further we have to specify the state \S) of the system. We take all 
electron and positron scattering states to be empty and the bound state is also unfilled 
(in other words this state coincides with the vacuum state as defined in ( |3.4| ). Then 

a e , (k) \S) = c e . {k) \S) = (3.12) 

b\S) = (3.13) 
According to ( |3.12| ) (S\ Q norm | -S) = and so 

(S\Q\S) =Q ,tmt (3.14) 

We now slowly increase the potential strength until eventually the bound state energy 
is found belowi? = 0. We denote everything pertaining to this later time by a ("). 
According to ( |3.1| ) the field now evolves as 

ij>(x, t) = J2k{ae(k, 0)«(+) e (A;, x)e~ iEt + a (k, 0)u (+)o (k, x)e~ iEt + 

(3.15) 

+cl(k, 0)u He (k, x)e iEt + c\{k, 0)u ( _ )o (A;, x)e iEt } + S (0)u(x)e iEbt 

Since the change in potential strength is assumed to be sufficiently slow and small and 
since there was no electron occupying the bound state initiallly, it is clear that there 
will still be no electron occupying the bound state. Nevertheless, by the conventional 
interpretation of Dirac's theory, the absence of an electron in a state with E < 
signifies the presence of a positron which in this case is a bound positron. For a slow 
enough change in the potential strength it is also clear that no radiation is emitted. 
Thus in terms of the (") operators the state \S) satisfies 

ae, {k) \S) = c e>0 {k) \S) = 



cfld\S) = \S) (3.16) 

Despite the appearance of a bound positron total charge is conserved, i.e. (S\ Q \ S) 
is still equal to Qo t init when Q is expressed in terms of the (") operators. To see 
this note that the d)d term in Q norm ( |3.7|) now contributes a term +1 (the positron 
charge) because of ( |3.16p . This however is cancelled by a term —1 originating from 
Qo simply by counting states; since one bound state has now dived below zero the 
vacuum asymmetry is equal to Qo,init — 1- 

Finally it may be of interest to examine the expectation value of the total charge 
when, as the potential changes, the energy of the bound state exactly reaches zero. 
We can then use either expansion 

ip(x, t) = J2k{a e (k, 0)u {+)e (k, x)e~ iEt + a a (k, 0)u {+)o (k, x)e~ iEt + 

(3.17) 

+4(k, 0)u He (k, x)e lEt + cl(k, 0)u Ho (k, x)e iEt } + bu{x) 
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or 

i/j(x, t) = J2k{ae{k, 0)u {+)e (k, x)e~ iEt + a D (k, 0)u(+) o (k, x)e~ iEt + 

(3.18) 

+4(k, 0)u ( _ )e (fc, x)e lEt + cl(k, 0)u Ho (k, x)e lEt } + d)u{x) 

The use of b or gF in ( p.!7|) and (|3.18| ) respectively is purely a matter of taste since 
in the absence of an energy exponential our convention dictates nothing. Assuming 
again that the bound state level has an electron vacancy ( |3.17|) entails 



b\S) = (3.19) 



whereas (|3.18|) entails 

tfd\S) = \S) (3.20) 
In both cases scattering states are vacant 

a e , (k) \S) = c e , (k) \S) = 

Notice however that the two conventions have different implications for the value of Qo 
when the level is at exactly zero energy. According to ( |3.19| ) the zero mode is counted 

as an electron state hence Qo of ( |3.8| ) gets a +- contribution due to normal ordering 

of b and . In contrast according to ( |3.20| ) the zero mode is counted as a positron 

state hence Qo gets a — — contribution due to normal ordering of d and d'. However 

physical results do not change. When we calculate (S\ Qnorm \S) then according to 
the first convention there is no contribution from the zero mode because of ( p.19 ). On 



the other hand according to the second convention there is a +1 contribution because 
of (|3.20|) and the presence of the d^d term in Qnorm- Thus the total contribution to 

(S\ Q \ S) due to the zero mode is +- as before. Notice that when the bound state 

is slightly above zero then its contribution to Qq is +- (according to ( |3.8| ) and its 

contribution to (S\ Qnorm \S) is zero according to ( |3.13|) . Thus the contribution of the 
bound state to (S\Q\S) depends smoothly on the potential. 
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4. Phase Shifts at Threshold and Levinson's Theorem and the 
Vacuum Charge. 

We now discuss Levinson's theorem [13], [14] , [15], [16], for a one- dimensional system 
since it is intimately related to the discussion of the vacuum charge on which our in- 
terpretation of positron radiation is based. We follow the strategy used by Barton [12] 
in his discussion of Levinson's theorem for the Schroedinger equation in one dimension 
and enclose the system in a box of length 2L with periodic boundary conditions 

^(-L)=^(L) (4.1) 

(Note that if one compares what follows with reference 12 one should interchange the 
roles of integers n and v.) We consider even and odd states separately. An odd state 
of positive energy has asymptotically the form 

, . . ( ik sin (kx ± A G+ ) \ , . 

i>oM (x - ±00) = ^ (£ _ m) ^ {kx ± + A ; +) J (4.2) 

and an even one 

, , ( ik cos (kx ± A e+ ) \ , . 

i> even (x -+ ±oo) = (e _ m) s . n {kx ± Ag+) j (4.3) 

where we neglect normalization factors (the subscript ± denotes the energy sign). 
Negative energy states are obtained by replacing e by — e as already mentioned in 
section 1. The phase shifts introduced in section 2 are connected to the above by 

5 + (e) = A e+ (e) + A 0+ (e) (4.4) 



5_ (e) = A e _ (e) + A _ (e) (4.5) 

Condition (4.1) implies 

k (i/ e ,±) L + A e± (k u ) = z/ e , ± 7r (4.6) 

k (i/ 0)± ) L + A Q± (fc„) = z/ ,±7T (4.7) 

v being integer. It is clear from (4.2), (4.3) that the wavevectors k are strictly positive; 
this will impose restrictions on v. In the continuum limit calculations similar to the 
ones outlined in section 2 yield 

A e (±m) = | + n e , ± (0)7r (4.8) 
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A (±m) = n ,± (0) 7T (4.9) 

(the argument in the above relations indictes that V — > 0). In the free case relations 
(4.6), (4.7) are still valid with no A appearing. For V small enough, the system reduces 
to that of the Schroedinger equation discussed in [12]. In that case Barton shows that 
for counting purposes it is sufficient to consider the number of positive energy states 
for a free particle to be integers v e ,+, ^o,+ ranging from a minimum value 

^minjree = 1 (4.10) 

to some large cut-off value N. We repeat the process for negative energy states. We 
thus have 2N even solutions and 2N odd ones (counting both positive and negative 
energy states). Suppose that we now switch on the potential keeping it arbitrarily 
small. We know that the number of positive energy even scattering states decreases 
by 1 since there is one even bound state appearing. We have the same number of odd 
scattering states and the same number of negative energy scattering states. The total 
number of states is conserved. In other words 

^min,o,— ^min,o,+ ^min,e,— 1 (4-11) 



^min,e,+ = 2 (4.12) 

It should also be observed that for any value of z/ m i n the wavevector k (f m m) goes 
to zero in the large L limit. We are thus entitled to replace A e o (k (iVin,±)) by the 
appropriate threshold value A (±m). Substituting (4.8) in (4.6) we get 

7T 

k (i/ e ,+) L = — — — n e>+ (0)7r + z/ e>+ 7r (4.13) 

k (i/ e _) L = ~- n e ,_ (0) 7t + z/ e ,_7r (4.14) 

Similarly 

k iyo,+) L = —n ^ + (0) 7r + f 0)+ 7r (4.15) 

k {v -) L — —n - (0) 7r + Vo-'K (4.16) 
From (4.11), (4.12) and from the fact that k is strictly positive we deduce that 

n e , + (0) = l (4.17) 

n e ,_(0) = (4.18) 
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In other words 



(m) = — , A e (-m) 



7T 



2 



Via a similar argument 



n ,+ (0) 



n ,_ (0) = 



and 



A (m) = A c (-m) = 



(4.19) 



This concludes the question of the determination of phase shifts at threshold mentioned 

7T 

after (4.5). Note that A e (±m) has a discontinuity of — at V = just as it does for the 

one dimensional Schroedinger equation since an attractive one dimensional potential 
no matter how weak always has at least one bound state. 

We can now turn on the potential gradually: then for specific values of V (cf. 
the discussion at the end of section 2) a scattering state (alternatively even or odd) 
crosses m and becomes bound. Since one scattering state is lost it is clear that v m i n ,e,+ 
or f m in,o,+ respectively increases by 1. From the requirement that k be positive we 
deduce from (4.11) and (4.17) that n ej+ (V) and n G>+ (V) denote the number of even 
and odd states respectively that have crossed m for the given value of the potential 
(including the even bound state that exists just below m for an arbitrarily small value 
of the potential). Similarly as the potential deepens further, bound states cross — m 
and join the continuum. The fact that more negative energy scattering states become 
available means that iVi n , e ,- anc ^ u mm,o- decrease and from (4.18), (4.16) we deduce 
that n 6i _ and n Gi _ decrease. Since (cf. (4.16), (4.18)) they start at zero for vanishing 
potential we conclude that n 6i _ and n 0j _ are negative and that the absolute values 
\n e - (V)\ and \n - (V)\ represent the number of even and odd (ex) bound states that 
have crossed — m for the particular value of V. Simple bookkeeping then yields that 
n et+ (V) + n e - {V) and n Qi+ (V) + n 0j _ (V) denote the number of even and odd bound 
states respectively. From (4.8) and (4.9) we get 



The above relations constitute Levinson's theorem in the present problem. 

The connection between the number of bound states and the jumps of phase shifts 
at threshold by n can be seen directly in the case of the square well examined in 
section 2. Return to expression (2.19) for the phase shift, set k ~ 0, E ~ m so that 

the quantity — in (2.19) tends to infinity regardless of the strength of the potential 

7 

and vary V. For some particular value V new a new bound state crosses E = m. Recall 
that according to the analysis of section 2.4 tan 2pa as a function of V vanishes at 




(A e (m, V) + A e (-m,V))-l 



(4.20) 



TI>bound,odd 



-(A (m,V) + A (-m,V)) 



(4.21) 



16 



V new . From (2.7) it is also clear that tan 2pa is an increasing function of V, so since 
V is increasing, tan 2pa is crosses zero from below. Thus according to (2.19) 

S (m, V new — 0) = arctan (— oo) 



5 (m, V new + 0) = arctan (oo) 

hence when a new bound state appears 5 (m) jumps by ir, i.e. n in (2.22) increases 
by 1. In other words n counts the number of bound states that appear (including 
the bound state that exists just below m for an arbitrarily small V). By a similar 
reasoning when a bound state crosses zero and disappears (for increasing V), 5 (—m) 
jumps by — tt, i.e. n! decreases by 1. Hence for an attractive potential n' is negative 
and its absolute value counts the number of bound states that disappeared. 

Before we move to the determination of the vacuum charge let us point out the 
constraints that (4.7), (4.9) impose on the minimum value of the wavevector k (for L 
large but finite). Setting A = in the above relations we get that for a free field in 
all cases (positive or negative energy, even or odd parity) 

TT 

km'm,free = ~j~ (4.22) 

corresponding to v — 1. For a very weak potential (4.7), (4.9) in conjunction with 
(4.15), (4.16), (4.19) yield (again in all cases) 

fcmin = ^ (4.23) 

To obtain the expression for the vacuum charge we convert sums over states using the 
standard expression that follows from (4.2), (4.3) 

1 r°° ( dS\ 
^(states) = =-J dkU + ^A (4.24) 

where the lower limit will be stated presently. Separating the contributions of odd 
and even states (3.8) gives 



^ 1 f 1 f°° It d^e+(E)\ 

Qo — dk I L H — + nf, ottnc ; even _| 

2 n Jk m s„ \ dk 




1 f°° ( d5 0+ (—E)\ If 00 ( 

+- / dk[L-\ °— + n houndMd& / dk \L + 

7T Jk min \ dk ) TT Jk min \ 

I f° Ale ( T dA °-( e "> \ _ 

^bound,even,— I dk I Li -\- I fl\) 0und odd,— 

TT Jk min \ dk J 

with /c m i n given by (4.23). It should be noticed that although k m i n tends to zero in the 
L — > oo limit this is counterbalanced in the integral by the factor L. Although this is 
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irrelevant for our purposes it is in general essential if we want to ensure conservation 
of the total number of states. The fact that the lower limits of integration in (4.25) 
are all identical allows the L proportional terms to cancel out after replacing oo by 
a cut-off wavevector K (this is true even after a more careful regularization. The 
dk integration over the A functions is performed trivially and brings in the functions 
S + (e), 5 -(e) via (4.4), (4.5). For convenience we define the number of bound states 
according to the energy sign and regardless of parity 

= flbound,even,+ Abound, odd, + 



N— — fl>bound,even,— ~\~ fl>bound,odd,— 

Then Levinson's theorem can be trivially rewritten 

N + + AL = - (S+(m) + <L(m)) - 1 

7T 

We finally get 

Qo = \{l (<M<x>, V) ~ 5 + K V) - <L(oo, V) + <L(m, V)) + N + - AL J (4.26) 
In the particular case of a square well (2.21), (4.26) yield 

Qo = U— + ~ (~h (m, V) + 5_ (m, V)) + N + -N\ (4.27) 

Z l 71 71 ) 

It may be appropriate at this point to enlarge on the role played by k m \ n introducted 
in (4.22), (4.23) and explain the meaning of the statement following (4.25). Suppose 
that the potential is extremely weak and that we have one (just) bound state. It is 
then natural to assume that the total number of positive energy states (scattering plus 
bound) is the same as before. Equation (4.24) gives 

L r K 1 r K dS 
(positive states, V -> 0) = - / dk + - dk— + 1 = (4.28) 

TT Jfc m ; n 7T Jk m i n UiK 

Lk L 1 

^min + {k m in) + 1 
71 71 71 

For a sufficiently weak potential 5 (K) can be neglected. As already explained for 

37T 

large LS (k min ) ~ S (0) — — . If we were to neglect the role of k min (i.e. set k m - m = 0) 
we would get 

LK 1 

^(positive states, V -> 0) = (4.29) 
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Given the first term in (4.29) would be there even in the free case it is clear that the 
above result and the idea of dropping /c min are absurd: (4.29) implies that switching 
on the potential creates half a state. However if we take (4.23) into account we get 

ELK L tt 1 3tt LK 
(positive states, V -> 0) = . — .— + 1 = 1 (4.30) 

TT TT 2L TT 2 TT 

On the other hand for strictly zero potential strength k m i n is given by (4.22), hence 

LK L tt TT LK 



y^Xpositive states, V = 0) 



TT TT L L TT 



in manifest agreement with (4.30). 

Let us focus on the behaviour of Qq for small changes of the potential. When the 
potential changes slightly and a new bound state appears the phase shifts <5(±oo) = 
±2Va change smoothly. The jump of 5(m) by tt in ( |4.26| ) is counterbalanced by the 
increase in N + by 1. Thus Qq also changes smoothly. On the other hand if for a 
slight change of the potential a bound state crosses zero then the phase shifts behave 
smoothly, N + decreases by 1 and AL increases by 1. Hence Qq decreases abruptly by 
1. This however should not necessarily be construed as a physical discontinuity. Recall 
for example the argument after ( |3.16| ): For the particular state \S) considered there 
the decrease in Qq is counterbalanced by a corresponding increase in (S\ Qnorm \ S) 
since the empty bound electron state now appears as an occupied positron state with 
charge +1. Thus the total charge has no discontinuity. 

Let us see this last point in somewhat greater detail. Suppose that one discrete 
level (the lowest lying even state in the present case) crosses zero for V = Vq. (For 
definiteness assume that a is such that the first odd state has already appeared; it 
can be checked that this is indeed feasible for say a = .7.) Then at V — Vq— the even 
bound state is just above zero and we assume it to be empty. The vacuum charge is 
obtained from (4.27) with N + = 2, iV_ = 0, n = 1, n' = 

Qo = ^ (4.31) 

TT 

and the expectation value of the total charge is 

(S\ Q(Vq-) \S) = (S\ Q norm (V -) \S) + Qq 

where the argument (Vq— ) serves to remind that creation and annihilation operators 
depend on V. Consider a very small change in the potential so that V = V +. Since 
the even bound state has dived below zero it is now a positron state and according to 
the hole theory (and to our conventions) the absence of an electron amounts to the 
presence of a positron. So in terms of operators pertaining to V — Vq+ the state \S) 
defined by ( |3.12j ) satisfies 

d^V +)d(V +) \S) = \S) (4.32) 
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Since the potential changed infinitesimally no radiation has been emitted. Hence 
according to (|3.7|) 



(S\ QnorUV + ) \S) = (S\ QnormiV^) \S) + 1 

where the last term is due to the presence of d\ ven d even in Q norm and simply reflects 
the fact that the positron carries charge +1. The important point is that charge is 
conserved. Indeed now iV + = 1, iV_ = 1 and according to (??) Q has decreased by 1: 

~ 2Vba 

Qo=— ~ 1 (4.33) 

7T 

The above reflects the connection between vacuum charge and spectral asymmetry. 
The charge expectation value (S\Q\S) rests unchanged during the crossing. 
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5. An Example: The 5 function Potential. 

We consider the attractive (for electrons) potential —X5(x). We define this potential as 
a square well in the limit V — > oo, a —>■ so that Va is kept finite and A = IV a. (This 
definition avoids the problems of definition discussed in reference [14]. The advantage 
of the 5 function limit lies in the fact that we can obtain simpler closed formulae for 
some of the quantities treated previously. 

We first turn to the determination of the scattering amplitude and the phase shift. 
For a wave incident from the left 

*<*<°> = (*- m )*""+*( ^ m )«"~ f") 

i>(x > 0) = F ( E * ^e"" (5.2) 
A solution to ( |2.6|) is clearly furnished by 

ip(x) = / dxexp(iu2(E + X5(x) — m))ip(xo) (5.3) 

Jx 

for some arbitrary xq. Applying the above equation for xq = 0—,x = 0+ we connect 

(PD, Q 

^(0+) = e lXa2 ip(0-) (5.4) 

or 

( * \+R\ j7 lk ) = (cos A + i sin Acr 2 )F ( ^ ) 
y E-m I \E -m ) y ' \ E - m J 

We thus get two equations for F, R. After trivial algebra 

k 



F(E) 



k cos A — iE sin A 

This leads to 



' E 

6(E) = arctan ( — tan A ) (5.5) 



k 

Let us compare the above result to (|2.19p . To this end the V — > oo limit in section 
2 should be taken from the start. Then clearly 2pa = A. Again in the same limit and 
after some algebra we get 

1+7 2 ^ E 
27 k 

Thus ( 2.19|) and ( |5 . 5| ) agree. In the high energy limit clearly 

5 (±00) = ±A (5.6) 
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(note that this agrees with the convention ( |2.20| )). It also agrees with result ( |2.21| ). 
At threshold we get 



5{m) 



TT 



+ n(\)% 



(5.7) 



8(— m) 



TT 



2 + n '( X > 



(5.8) 
at 



The analysis of the previous section still holds and entails n(0) = n'(0) 
threshold. 

We now turn to the spectrum of the bound states. The wavefunctions are of the 
form 



— K 

E — m 



-k\x\ 



modulo a sign depending on parity. Applying fl5T4|) we get for even bound states 



— K 

E — m 



cos A sin A 
— sin A cos A 



— K 

E — m 



(5.9) 



and for odd bound states 

( - K 

\ E — m 



cos A sin A 
— sin A cos A 



— K 

E — m 



(5.10) 



Using relation fl2.10|) between k and E we get the equations determining the spectrum. 
For even bound states 

A m — E 



tan ■ 



and for odd bound states 



tan ■ 



A 



m + E 



lm + E 
m — E 



(5.11) 



(5.12) 



Equations fl5.11| ), ( |5.12|) are identical to ( |2.11| ), (??) in the V — > oo limit. Both give 
the equation for the spectrum 



E = m cos Asign(sin A) 



(5.13) 



A A 7T 

Thus even bound states appear when — = Ntt and disappear when — = — + Ntt. 

A 7T A 

Odd bound states appear when — = — — + Ntt and disappear when — = Ntt. A 

distinctive feature of the 8 function potential is that whenever a bound state of a 
certain parity appears a state of the opposite parity disappears (and vice versa). For 
example when A = tt the original (even) bound disappears and the first odd one 
appears, at A = 2tt the odd state disappears and the second even state appears etc.. 
In other words Amod7r counts the number of bound states (excluding the first one that 
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exists for arbitrarily small A) that have crossed E = m. These results may also be 
obtained from the analysis of section 2.4 in the limit a — > (thus essentially dropping 
m everywhere). Finally from ( |5.13| ) we deduce that a discrete level crosses zero when 

A = — + Ntt. 
2 

Let us examine the vacuum charge in the 6 function case. Suppose that iV bound 
states (on top of the original half bound state existing for arbitrarily small A) have 
crossed m, i.e. A = Ntt + A', 0< A' < tt. Then <5(oo) = A , S(— oo) = —A, 6(m) = 

— + Ntt, 8(—m) = — — — Nn. If the bound state is above zero, i.e. iV + = 1 , AL = 
flOp gives 

Qo = - (5.14) 

7T 



In the same way if N + = , N_ 



Qo = --l (5.15) 

7T 



TT 1 

Thus at A = — (when the crossing takes place) the vacuum carge jumps from +— to 
— , the difference being correctly equal to 1. As explained towards the end of section 

TT 

3.1 the value of the vacuum charge at A = — depends on our convention concerning 
the labelling of the zero mode. 
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6. The Transition from a Subcritical to a Supercritical Poten- 
tial and Positron Emission. 

We now finally consider the transition from a subcritical to a supercritical potential. 
We assume that the potential starts at a slightly subcritical value V su b, makes an 
abrupt change to a supercritical value V super at time t\, stays fixed at this value until 
time ti and then makes an abrupt jump back to V su b. By assumption there is a positron 
occupying the bound state for if we had started with an initially vacant electron bound 
state (for a weak potential) then as the potential became stronger the bound state 
would have crossed E — and then the absence of an electron in that state is now 
interpreted as the presence of a positron. The implicit assumption is that the initial 
electron vacancy (or equivalently positron prescence) persists during the switching on 
process from V wea k to V su i>. Wavefunctions pertaining to the supercritical potential 
and corresponding creation and annihilation operators will be denoted by a (~), the 
time argument of the latter ranging from t\ to £2- 

The crucial observation is that although the potential may change abruptly at 
t = ti the field ip(x, t) is continuous. Thus we can expand ip(x, t±) in terms of operators 
pertaining both to t — ti— and to t — t\ + . For t — ti— the expansions read 

ip(x, ti) = J2ki a e(k, ti)u(+) e (k, x) + a a (k, ii)u( +)o (A;, x)+ 

+4(k, ti)M ( _ )e (/c, x) + cl(k, ti)u(- )o (k, x)} + d\ ven {ti)u even (x) + & ocW (ii)u dd(£i) 

(6.1) 

ip^(x, h) = Eki a l( k , ti)u\ +)e {k, x) + al(k, h)u[ +)o (k, x)+ 

+c e (k, *i)«J_ )e (A;, x) + c (k, *i)«J_ )o (A;, x)} + d even {ti)u even {x) + &Ld(^i) n odd(^i) 

(6.2) 

and for t — t±+ 

ip(x, ti) = J2k{a e (k, ti)u(+) e (k, x) + a (k, ti)u^ +)o (k, x)+ 

(6.3) 

+4(k, i 1 )«(_ )e (A;, x) + cl(k, ti)M ( _ )o (/c, x)} + &odd(ii)M od d(ti) 
ip*(x, h) = E k {at(k, *i)u[ +)e (fc, x) + al(k, h)u\ +)o (k, x)+ 

(6.4) 

+c e (k, *i)«J_ )e (fc, x) + c (k, ii)«J_ )o (A;, x)} + 6Lj(*i)«Lj(*i) 
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Expansion ( |6.3|) can be inverted to give 



a e (k,ti) = I dxu^ e (k,x)ip(x,ti) (6.5) 



$ e (k,ti) = / dxu\_s (k,x)ip(x,t 1 ) (6.6) 



Similar expressions hold for the odd operators and for the bound state ones. Ex- 
pressions for the Hermitean conjugates can be obtained by taking the complex con- 
jugates of the right hand sides of (|6.5|,|6.6|). We invoke continuity and substitute (|6.3|) in 



75]), (|6.6|) . We notice that we get overlap integrals of the form / dxu\±y a (k, x)u[±) e>0 {k' 



x ) 



Given that outside the well the wavefuctions are plane waves such integrals are pro- 
portional to 5 kk >5 eo . We thus get 

a e (k, h) = A k a e (k, h) + B k cl(k, h) + F k d) even {t x ) (6.7) 
c\{k, h) = G k a e (k, h) + L k c\{k, t x ) + MkdlM (6.8) 

where 

A k = f dxul +)e (k,x)u(+) e (k,x),B k = J dxu\ +)e (k, x)u(-) e (k, x) 
F k = J dxu\ +)e (k, x)u even (x), G k = J dxu\_ )e (k, x)u( + ) e {k, x) (6.9) 

L k = Jdxu\_ )e (k,x)u He (k,x), M k = J dxu\_ )e (k, x)u even (x) 

Expressions for a\(k, t x ) and c e (k, tx) are obtained by taking the Hermitean conjugates 
of (|6.7p , ( |6.8|) respectively. The corresponding equations for the odd modes of course 
involve the odd rather than the even bound state because of parity 

a (k, h) = A' k a (k, t x ) + B' k cl(k, h) + F' k b oM {t x ) (6.10) 

cl(k, t x ) = G' k a (k, t x ) + L' k cl(k, t x ) + M' k b odd (t x ) (6.11) 

where the coefficients A' k , B' k , F k ,G' k , L' k , M' k are defined as in ( |6.9|) replacing even by 
odd. For the anticommutation relations 

{a e (k,t),al(k',t)} = 5 kk > 
{c e (k,t),cl(k',t)} = 8 kk > 
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to be valid we should have 

A* k A k + B* k B k + F*F k = l 

G* k G k + L* k L k + M* k M k = 1 

(no summation over k). Both of them are satisfied due to the fact that the set of 
wavefunctions u^ e (k,x),u^ e (k,x),u even (x) form an orthonormal set. Similar rela- 
tions are valid for the primed coefficients. We now make some rather drastic approx- 
imations. Since the potential is assumed initially to be just subcritical the even state 
lies very near the negative energy continuum. Also the change in the potential is 
assumed to be small so that we end up with a slightly supercritical potential. Thus it 
is natural to assume that there is negligible overlap between the new positive energy 
wavefunctions on the one hand and the old negative energy wavefunctions and the 
bound state on the other. Hence in ( |6.9| ) A k = B k = F k = and 



a e (A;,ti) = a e (Mi) (6-12) 

a+(Mi) = 4(Mi) ( 6 - 13 ) 
By the same token we can take G k = 0. Thus 

c e (k, t x ) = M* k d h {t x ) + L* k c e (k, t x ) (6.14) 

ct{k, h) = Mkdlih) + L k ct{k, h) (6.15) 



Note that in the approximation we are working (neglecting the overlap integrals of the 
bound state with new states lying in the positive energy continuum), i.e. essentially 
assuming completeness of the set u^ e (k, x), and using the orthogonality property of 
the above set we get 

£M fc M* = l (6.16) 

k 

Following a similar reasoning we deduce that all odd operators referring to t = ti- 
ara equal to the corresponding odd operators referring to t = t\+. 

Having obtained the relations beteen operators we come to the definition of the 
state IS") of the system in terms of operators pertaining to t = ti — . Recall that since 
we work in the Heisenberg picture the state vector is time independent however the 
vacuum (as all basis kets) is time dependent. The statement that there is one positron 
occupying the even state is equivalent to 

^(tO^ti) \S) = \S) (6.17) 

We also assume that there are no free electrons or positrons, i.e. 

a e , Q (k, h) \S) = c e , (k, ti) \S) = (6.18) 
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As already mentioned we can arrange the parameters so that when the first even state 
merges with the negative continuum only one (odd) bound state exists. Although it 
does not matter in what follows and simply in order to fix the notation we have to 
decide whether this odd state is above or below zero when the potential assumes its 
critical value. We take the first possibility; this is indeed the case for the value a = .7 
mentioned in the previous section. We assume that there is an electron vacancy in 
the odd state, i.e. 

bodd(h) \S) = (6.19) 
Relations Q6.17| ), ( |3.12| ), ( |3.13| ) specify the state completely. Using the above equa- 



tions together with (|6. 14 ) , ( 6.15Q we can calculate the average number of positrons 



with momentum k for the supercritical potential: 

N eM = (S\~ct(k)~c e (k)\S) = M k M* (6.20) 

(|6.16|) gives 



J2 N e,k = l (6.21) 



i.e. there is one positron emitted. Expression (|6.20 ) for N e ^ yields information on the 



energy spectrum of the positron emitted. Notice that the mixing between creation 
and annihilation operators at t — tx, the ordering of the operators in ( |6.20| ) and the 
use of ( |6.16i ) are crucial. 



It should be clear from the above discussion that if the original bound state is filled 
by an electron then 

d(h) \S) = 

and the above reasoning leads to the conclusion that no positron is emitted. 

We should check whether positron emission as described above is consistent with 
charge conservation. The Hamiltonian appropriate for the time dependent well consid- 
ered here is given by (|3.10|) hence the total charge defined by (|3.5|) is indeed conserved. 
To demonstrate this we have to evaluate the charge expectation value (S\ Q \S) with 
Q expressed in the form (|3.6|) in terms of the (~) operators pertaining to t±+ and 
compare it to (S\Q \S) evaluated at t = t\— . At t%+ 

Q = *i) + c\(k, t x )c e {k, *i) - al(k, *i)5 (A;, * x )+ 

(6.22) 

0,sup er 

and at t\— 

Q = Efc{-a e (^, h)a e (k, t t ) + c\{k, t^c^k, ti) - a\(k, h)a (k, t x )+ 

(6.23) 

+4(Ml)Co(Ml) - 6Ld(*l) & odd(*l) + dlveni^devenih)} + Qo,sub 
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Using (|6.21 ) we notice that the increase of the (?c contribution in the evaluation of the 



charge expectation value by one in ( |6.22j ) exactly balances the contribution of the d)d 
in fl6.23| ), thus the contribution of the normal ordered charge Q nor m rests unchanged. 
We also observe that Qo tSU b = Qo, S u P er when the potential becomes supercritical: An 
energy level crosses — m but this does not change the spectral asymmetry. (The latter 
would change only in the case where the energy level crosses zero but this does not 
happen at t = ti.) One can see this either from (|3.8|) by inspection or in more detail 
from (|4.26|) : when a bound state merges with the negative energy continuum AL 
decreases by 1 but S(— m) increases by ir, the net result being zero. 

Regardless of the experimental feasibility we can in principle arrange for the poten- 
tial to achieve the supercritical value and then remain fixed. We would then observe 
one positron emitted while V goes from V su h to V^ up er and nothing else. Thus the tran- 
sition from a subcritical to a supercritical potential creates just one positron in this 
case and not an electron-positron pair as Zeldovich and Popov (and others) predict. 
We shall see presently how and in what sense an additional electron appears. 

At time ti+ the positron charge density is concentrated inside the well and subse- 
quently leaks away. Since the ex bound state has been transformed to a transmission 
resonance the characteristic time scale for the process is related to the corresponding 
time delay. Since we intend to switch back to a subcritical potential at time t 2 the 
time interval t 2 — £1 should be rather large compared to the time delay so that the 
positron has time to escape. An estimate of the lifetime is provided in Appendix A. 

When we switch back to the subcritical value the positron has escaped from 
the well and its charge is carried by travelling waves, hence the expectation value 
(S\cl(k)c e (k)\S^ has increased by 1 compared to its value for t < t\. We have already 
shown that Qq stays unchanged when a level crosses the threshold —m. Thus charge 
conservation requires that there is no positron occupying the even bound state for 
t > t 2 - When the potential returns to its original value V su b the absence of a bound 
positron signifies the presence of a bound electron. Thus the round trip V su b — > K, uper 
— ► Vsub starting with an empty electron bound state results to the emission of one 
positron and the creation of one bound electron. It is only in that sense that one may 
talk about pair creation. It should be clear that the roles of positron and electron are 
rather unequal, the former being free and the latter being bound. Their roles would 
of course be reversed for a potential of a different sign. It must be stressed again that 
if the bound state were filled by a positron when the potential is slightly subcritical 
then nothing would happen when the potential becomes supercritical. 

If one prepares the next (higher lying) odd bound state to be empty then as the 
potential becomes even stronger, this state merges with the negative energy continuum 
when 

V = Vo = \ + m 2 + m 
V a 2 

and there is one further positron emitted. On the return trip the odd bound state is 
occupied by an electron. 
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7. Conclusions 



We hope that we have now shown how to write down a second quantised theory of 
spontaneous fermion production in the context of a Dirac particle bound by an ex- 
ternal field. Dirac's own understanding of a physical positron state as an unfilled 
electron state of negative energy is robust and perfectly capable of dealing with this 
phenomenon. Considerations of the vacuum charge (and the related topic of Levinson's 
theorem) are required for a full description of the problem. Zeldovich and Popov's 
intuition which allowed them to give a detailed analysis of spontaneous positron pro- 
duction based solely on the first quantised Dirac equation is generally sound. The one 
place it fails is in its identification of the cause of positron production with electron- 
positron pair production. This identification can be made when the supercritical 
potential is switched on at early times and off at late times, but when the potential 
is switched on and remains on then spontaneous positron production is distinct from 
electron-positron pair production. 

We have not shown why narrow positron peaks are observed experimentally in 
heavy ion collisions [17]; in principle the overlap integral (6.20) taken between re- 
alistic Coulomb wave functions should provide a clue. A detailed knowledge of the 
time-dependence of the external potential is, however, probably more important. In 
this paper, we only assume a trivial time dependence in order to demonstrate the 
fundamental principles of spontaneous positron radiation, but we hope to return to 
the more general problem of treating electron-positron pair production in a time- 
dependent external field in a forthcoming paper. 
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A. Appendix A. 



At 



According to nonrelativistic quantum mechanics the time delay At of a wavepacket 
built around a transmission resonance is given by 

1 d5 

vq dk 

where vq is the group velocity of the wavepacket, k the wavevector corresponding to 
the resonance energy and 5 the phase shift (see e.g. [10]). In the present case Vq is the 
relativistic group velocity 

k 



V 



\E\ 



E l -k l 



m 



The phase shift is given by expression (|2.19|) 



/1 + 7 2 \ 
6(E) = arctan I — tan2pal — 2ka 



(a.i; 



where 



k E + V + m 



7 = — 



p E + m 

and p = J[E~+ V) 2 — m 2 . We are interested in the N — 1 transmission resonance, 
hence p satisfies 

2pa = 7T (A.2) 



From (|A.1| ) we get 



d5_ 

dk 



a 7 + 



l\ dp 



2a 



(A.3) 



7/ dk 

where we have used ( |A.2|) . Using the expressions for p and k in terms of E we readily 
calculate the derivative 

dp _ k E + V 
dk p E 

Then after some algebra 



l\ dp E + V ( E — m 



^ ^ 7 J dk 



E 



E + m 



E +V -m E +V +m 



(A.4) 



The above expression is exact. We now make the approximation of a just supercritical 
potential. Hence we can replace E by — m and V by Vf as given by ( |2.31| ). Substituting 
([ATI ) in (|A~!) we deduce 

dS m 

-rr = 2a , 



dk 



4^2 



+ m 2 — m 
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